Do interactions increase or reduce the conductance of disordered electrons? 

It depends! 
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We investigate the influence of electron-electron interactions on the conductance of two-dimensional 
disordered spinless electrons. We present an efficient numerical method based on diagonalization in 
a truncated basis of Hartree-Fock states to determine with high accuracy the low-energy properties 
in the entire parameter space. We find that weak interactions increase the d.c. conductance in 
the strongly localized regime while they decrease the d.c. conductance for weak disorder. Strong 
interactions always decrease the conductance. We also study the localization of single-particle 
excitations at the Fermi energy which turns out to be only weakly influenced by the interactions. 
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The influence of electron-electron interactions on the 
transport in disordered electronic systems has been inves- 
tigated intensively within the last two decades . Re- 
cently, the problem has reattracted a lot of attention af- 
ter experimental |^ and theoretical Q results challenged 
established opinions. 

It is well accepted ^ that non-interacting elec- 
trons in three dimensions (3D) undergo a localization- 
delocalization transition at finite disorder. In contrast, 
all states are localized in 2D and ID even for infinites- 
imal weak disorder However, today it is believed 
that the metal-insulator transition (MIT) in most ex- 
perimental systems cannot be explained based on non- 
interacting electrons. The metallic phase of disordered 
interacting electrons has been studied intensively within 
the perturbative renormalization group (RG) Q leading 
to a qualitative analysis of the MIT and the identifica- 
tion of different universality classes. One of the results 
is that the lower critical dimension of the MIT is c?^ = 2 
as it is for non-interacting electrons. Therefore it came 
as a surprise when experiments |^ on Si-MOSFETs re- 
vealed indications of a MIT in 2D. Since these experi- 
ments are performed at low electron density where the 
Coulomb interaction is particularly strong compared to 
the Fermi energy, interaction effects are a likely reason for 
this MIT. A complete understanding has, however, not 
yet been obtained. Explanations were suggested based 
on the perturbative RG @], non-perturbative effects d), 
or the transition being a superconductor-insulator tran- 
sition rather than a MIT [§. 

Theoretically, surprising results have been obtained for 
just two interacting particles in the insulating regime |Q. 
It was found that two particles can form a pair whose 
localization length is much larger than that of a single 
particle. Later an even larger delocalization was sug- 
gested for clusters of three or more particles ||T^. In 
the case of a repulsive electron-electron these delocalized 
states have rather high energy, thus their relevance for 
the low-energy properties of a degenerate system is not 
clear. It has been argued that the many-particle problem 



FIG. 1. d.c. conductance G(0) for a system of 5^ sites as a 
function of interaction strength U for different kinetic energies 
t. The disorder is fixed &tW=l. The statistical accuracy is 
better than the symbol size. 



can be reduced to a few interacting quasiparticles above 
the Fermi surface |ll| ]. This is, however, only possible, if 
the interactions do not change the nature of the ground 
state. All in all, not even the qualitative influence of 
interactions is understood in the insulating regime. 

We have numerically studied disordered 2D spinless 
electrons. Our calculations are summarized in Fig. ^ 
which is the main result of this Letter. It shows that the 
influence of repulsive electron-electron interactions on the 
d.c. conductance is opposite for high and low kinetic en- 
ergies (i.e. weak vs. strong disorder). The conductance 
of strongly localized samples {t — 0.01 to 0.03) is con- 
siderably enhanced by a weak interaction. With increas- 
ing kinetic energy the relative enhancement decreases as 
does the interaction range where the enhancement oc- 
curs. The conductance of samples with the highest ki- 
netic energies {t = 0.3 and 0.5) is reduced even by weak 
interactions. In contrast, sufficiently strong interactions 
always reduce the conductance. This is not surprising 
since for large enough interaction strength the system 
will form a Wigner glass. 

These findings shed some light on seemingly contra- 
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dieting numerical results on the transport of disordered 
spinless electrons in the literature. Studies (l2j of a 2D 
model in the diffusive regime yielded that interactions de- 
crease the conductance. The same conclusion was drawn 
from density-matrix RG studies and exact diagonal- 
izations in ID. In contrast, for 2D models in the 
localized regime p5| , pj[ | it was found that interactions 
lead to a delocalization. Up to now it has been un- 
clear whether these inconsistent results are due to be- 
ing in different parameter regions (weak vs. strong dis- 
order), different quantities studied (conductance, many- 
particle level statistics or charge stiffness), or long-range 
vs. short-range interactions. The results of this Letter 
suggest that being in different parameter regions is the 
most likely reason for the differences between the results 
cited above. A result similar to ours was obtained re- 
cently in a study of the ground state phase sensitiv- 
ity in ID. It was found that for small disorder repulsive 
(nearest-neighbor) interactions reduce the phase sensi- 
tivity while for large disorder the phase sensitivity shows 
pronounced peaks at certain values of the interaction. 

In the remainder of the Letter we explain the model 
and the calculational method and further discuss the re- 
sults. We consider a 2D quantum Coulomb glass model 
| p5| , p^ 18 1^. It is defined on a square lattice with 
M — sites occupied by = KM spinless electrons 
{0<K <1). To ensure charge neutrality each site carries 
a compensating charge of Ke. The Hamiltonian reads 



-I- lp^n^ - 



(1) 



where cj and Ci are the creation and annihilation oper- 
ators at site i, rii 
nearest neighbors. 
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and {ij) denotes all pairs of 
/vij represents the Coulomb 
interaction which is parametrized by its nearest-neighbor 
value U and t is the kinetic energy. The random poten- 
tial values Lpi are chosen from a box distribution of width 
2W and zero mean. (We always set W — 1.) Two im- 
portant limiting cases of the quantum Coulomb glass are 
the Anderson model of localization (for [/ = 0) and the 
classical Coulomb glass (for t = 0). 

The numerical simulation of disordered many-particle 
systems is one of the most complicated problems in com- 
putational physics. First, the size of the Hilbert space 
grows exponentially with the system size making exact 
diagonalizations of the Hamiltonian impossible already 
for very small systems. Second, the presence of disor- 
der requires the simulation of many different samples to 
obtain averages or distributions of physical quantities. 
For disordered interacting electrons the problem is made 
worse by the long range of the Coulomb interaction which 
has to be retained for a correct description of the insu- 
lating phase where screening breaks down. 



In this Letter we suggest an efficient numerical method 
to simulate disordered interacting electrons. It is based 
on the idea of configuration interaction adapted to 
disordered lattice models. The method, which we call the 
Hartree-Fock based diagonalization (HFD), consists of 3 
steps: (i) solve the Hartree-Fock (HF) approximation of 
the Hamiltonian as in Ref. ||l^, (ii) use a Monte-Carlo 
algorithm to find the low-energy many-particle HF states 
(Slater determinants), (iii) diagonalize the Hamiltonian 
in the basis formed by these states The HF basis 

states are comparatively close in character to the exact 
eigenstates in the entire parameter space. Thus it is suffi- 
cient to keep only a small fraction of the Hilbert space to 
obtain low-energy quantities with an accuracy compara- 
ble to that of exact diagonalization. The HFD method is 
very flexible, it works well in any spatial dimension, and 
is capable of handling long-range and short-range inter- 
actions. A detailed description will be given elsewhere. 
Most of our calculations have been performed for lattices 
with 5^ sites and 12 electrons keeping 500 basis states. 
We used periodic boundary conditions and the minimum 
image convention. We also studied 4^ and 6^ systems 
with K = 0.25 and 0.5 keeping up to 2000 out of 9 * 10^ 
basis states. 

We now turn to the conductance which we compute 
from linear response theory. The real (dissipative) part of 
the conductance (in units of e^ /h) is given by the Kubo- 
Greenwood formula |22|, 



Re G^-(c^) = ^ y^\{Q\f\v)\H{u + Eo- E,) (2) 

where lo denotes the frequency, is the x component 
of the current operator and v denotes the eigenstates of 
the Hamiltonian. Eq. (^) describes an isolated system 
while in a real d.c. transport experiment the sample is 
connected to contacts and leads. This results in a finite 
life time r of the eigenstates leading to an inhomoge- 
neous broadening 7 = r^^ of the 8 functions in (|2|) P3| . 
To suppress the discreteness of the spectrum of a finite 
system, 7 should be at least of the order of the single- 
particle level spacing. For our systems this requires a 
comparatively large 7 > 0.05. We tested different 7 and 
found that the conductance values depend on 7 but the 
qualitative results do not pi] . 

In a random system different samples will have differ- 
ent conductance values. Fig. |^ shows the probability dis- 
tribution P[log(G(0))] for systems in the localized regime 
with and without interactions. Both distributions show 
the same qualitative behavior, they are close to normal 
distributions corresponding to very broad distributions 
of the conductances themselves. The arithmetic average 
of the conductance is therefore not a good measure of the 
typical behavior. We instead use the logarithmic (i.e. ge- 
ometrical) average Gtyp — exp(log(G)) |^^, usually over 
400 disorder configurations. 



2 




-4 -2 

log (5(0)) 

FIG. 2. P[log(G(0)] for W = 1, t = 0.1, and 7 = 0.05. The 
histograms represent 2000 samples. The smooth lines are fits 
to Gaussians. The data for U — 1 have been shifted by 0.2 
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FIG. 3. G{lu) iorW = l,t = 0.03, 7 = 0.05. 

In Figs. H and || we present results on the dependence 
of the conductance on the interaction for two sets of pa- 
rameters. In Fig. ^ the kinetic energy is very smaU 
{t = 0.03). Thus the system is in the strongly localized 
regime, as we also estimated from the single-particle par- 
ticipation number Pgp ~ 2. Here a weak Coulomb inter- 
action (U = 0.5) leads to an increase of the conductance 
at low frequencies. If the interaction becomes stronger 
the conductance decreases and finally {U = 2) falls be- 
low the value of non-interacting electrons. We emphasize 
that the increase of the conductance for weak interactions 
is a true correlation effect: Within the HF approximation 
p8[ interactions always lead to a decrease of the conduc- 
tance. The behavior is qualitatively different at higher 
kinetic energy {t = 0.3) as shown in Fig. ^. Here the 
system is approaching the diffusive regime (Pgp > 10). 
Already a weak interaction {U = 0.5) leads to a reduc- 
tion of the low-frequency conductance compared to non- 
interacting electrons. If the interaction becomes stronger 
the conductance is decreased further. We have performed 
analogous calculations for kinetic energies t = 0.01. ..0.5 
and interaction strengths U — 0...2. The resulting d.c. 
conductances are those presented in Fig. |l|. 



• 


u=o 




U=0.5 


V 


U=1 


□ 


U=2 




0.5 1 1.5 

CO 

FIG. 4. Same as Fig. | but for t = 0.3. 




FIG. 5. Comparison of G(0) for W = l,t 
ferent system sizes and filhng factors. 
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We also checked for system size and filling factor de- 
pendences by simulating systems with 4^ and 6^ sites and 
filling factor K — 0.25 in addition to 0.5. We found the 
qualitative picture (as presented in Fig. [^) to be the same 
in all cases. As an example. Fig. || shows the interaction 
dependence of G{0) for t — 0.01 for the different systems 
studied. Clearly, the interaction induced enhancement of 
the conductance exists in all cases. Moreover the relative 
enhancement seems to increase from the 4^ system to the 
6^ system. (A comparison of even and odd linear system 
sizes is problematic since at half filling a regular array of 
charges is impossible for odd sizes. Moreover, any quanti- 
tative comparison of different sizes would require a more 
realistic description of the broadening.) 

In order to find out to what extent the behavior of 
the conductance is reflected in single-particle localization 
properties we also computed the single-particle return 
probability 
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G,,{e^i5)G,j{E-i5) . (3) 



Here Gij{e) is the single-particle Greens function. Rp{e) 
is the generalization of the inverse participation num- 
ber Pgp^(e) (of a single-electron state) to a many-particle 
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FIG. 6. Rp{e) for W = l,t = 0.1. The data are averaged 
over 2000 disorder configurations (10000 for non-interacting 
electrons) . 

system. Fig. ^ shows a typical result for Rp{e). We 
performed analogous calculations for t = 0.01. ..0.5 and 
U = 0...2. For all cases we obtain the same qualitative 
behavior: Close to the Fermi energy the return proba- 
bility is only weakly influenced by the interaction. Di- 
rectly at the Fermi energy, which is not accessible in our 
simulations because of our still too small system sizes, 
there may develop a slight enhancement of the return 
probability as a result of the Coulomb gap in the single- 
particle density of states. Such an enhancement has al- 
ready been observed within the HF approximation 
Within the results obtained in this Letter, the effect, if 
any, is weaker than within HF. For energies away from 
the Fermi energy the single-particle excitations in the in- 
teracting system become strongly delocalized compared 
to the non-interacting case. The interaction dependence 
of the conductance discussed above is, however, not re- 
flected in the single-particle return probability. 

To summarize, we have used the Hartree-Fock based 
diagonalization (HFD) method to investigate the trans- 
port properties of disordered interacting spinless elec- 
trons. We have found that a weak Coulomb interaction 
can enhance the conductivity of localized samples con- 
siderably while it reduces the conductance in the case of 
weaker disorder. If the interaction becomes stronger it 
eventually reduces the conductance also in the localized 
regime. Let us finally mention that although we show 
that interactions can enhance the conductivity in certain 
parameter regions this does not directly provide an ex- 
planation for the MIT in 2D ||^ since the importance of 
the spin degrees of freedom for this transition is estab- 
lished experimentally |2^. We emphasize however, that 
our method is very easy to generalize to electrons with 
spin. Work in this direction is in progress. 

We acknowledge financial support by the Deutsche 
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